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The new phase of a gauge theory in which the instantons 
are "polarized" , i.e. have the preferred orientation is discussed. 
A class of gauge theories with the specific condensates of the 
scalar fields is considered. In these models there exists an in- 
teraction between instantons resulting from one-fermion loop 
correction. The interaction makes the identical orientation of 
instantons to be the most probable, permitting one to expect 
the system to undergo the phase transition into the state with 
polarized instantons. The existence of this phase is confirmed 
in the mean-field approximation in which there is the first or- 
der phase transition separating the "polarized phase" from the 
usual non-polarized one. The considered phase can be impor- 
tant for the description of gravity in the framework of the gauge 
field theory. 

PACS: 04.60, 12.25 



I. INTRODUCTION 

In this paper the new nontrivial vacuum state of gauge 
theory is discussed. In this state instantons are "polar- 
ized" , i.e. have the preferred orientation. This polar- 
ized state is interesting by itself, but its importance is 
enhanced by Ref. Q where a scenario was suggested to 
describe effects of gravity in the framework of Yang-Mills 
gauge theory. In this approach space-time is supposed to 
be basically flat, there are only usual for the Yang-Mills 
theory fields of spin 0, 1/2 and 1 on the basic level of the 
theory. There are neither gravitons nor the dimensional 
Newton gravitational constant in the basic Lagrangian of 
the theory. In spite of its "conventional nature" the gauge 
theory as discussed in can provide a description of the 
effect of gravity. In recent decades it was supposed that 
quite new basic physical conceptions, such as supergrav- 
ity, strings and superstrings, for a review see are vital 
for quantum gravity. 

In Ref. a specific phenomenon based upon nontriv- 
ial topological excitations of the gauge field was consid- 
ered. The most clear topological excitation is the instan- 
ton [||. It has several degrees of freedom: its position, 
radius and orientation. Generally speaking an instanton 
may be oriented arbitrarily. It was assumed in ||l| that it 
is possible to construct such the Yang-Mills theory that 
the instantons in the vacuum state have a preferred di- 
rection of their orientation. That means that the prob- 
ability for any instanton to be oriented along certain di- 
rection is greater than along any other direction. This 
phenomenon can be called "a polarization" of instantons. 
For formulation of the problem discussed in this paper it 
is important to recall more precisely the main result of . 
With this purpose consider 5*0(4) gauge theory. Then the 
gauge algebra so(4) is the sum of two su(2) subalgebras: 
so(4) = su{2) + su{2). It was assumed in that the vac- 
uum state of 50(4) gauge theory possesses the following 
properties. First, instantons belonging to one su{2) sub- 
algebra of the gauge algebra are polarized. At the same 



time the antiinstantons belonging to this su{2) subalgebra 
are not polarized. Second, there is the reversed situation 
in the other su{2) subalgebra: the antiinstantons belong- 
ing to that subalgebra are polarized while the instantons 
of that subalgebra are not polarized. In other words, in 
both su{2) subalgebras the situation is nontrivial: in one 
of them the instantons are polarized while in the other one 
the antiinstantons are polarized. Note however, that the 
topological charge in every one su{2) subalgebra might be 
zero resulting in equal concentrations of instantons and 
antiinstantons for the given subalgebra. 

The concentrations of these polarized instantons of one 
su{2) subalgebra and polarized antiinstantons of the other 
su{2) subalgebra are supposed to be finite and equal. One 
can call this phenomena the "condensate of polarized in- 
stantons and antiinstantons" or "Instanton- Antiinstanton 
Polarization" in the vacuum. Certainly this desired vac- 
uum state can not be achieved in the framework of the 
pure gauge theory. The theory must acquire some special 
properties in the scalar and fermion sectors as discussed 
in detail below. The main result of Ref. [|l| is that if the 
vacuum of the 50(4) gauge theory possesses the conden- 
sate discussed then this gauge theory reveals the effect of 
gravity. 

Note first of all that existence of the considered nontriv- 
ial phase of the gauge field does not come into contradic- 
tion with the gauge invariance of the theory. The gauge 
invariance forbids those nontrivial phases whose order pa- 
rameter is not invariant under local gauge transformations 
Q , 1^ , 1^ . The orientation of instantons cannot be var- 
ied by the local gauge transformations. This is the well 
known fact for the pure gauge theory, see Refs. §, §. 
Therefore the instanton orientation may play the role of 
an order parameter for the considered nontrivial phase of 
the gauge theory. 

The most "natural" way to look for the condensation 
of polarized instantons is to find an interaction between 
instantons which forces any two instantons to have the 
identical orientation. Then one can expect the system to 
undergo a phase transition into the state with polarized 
instantons. An attempt to implement this idea meets two 
problems. First, in a pure gauge theory instantons do 
not interact. There is the well-known interaction between 
instantons and antiinstantons but the interaction be- 
tween instantons seems to be questionable because there 
is an exact multi-instanton solution with arbitrary orien- 
tation of instantons |ic| ]. Second, along with instanton 
interaction there must be no interaction between antiin- 
stantons in order to keep the system of antiinstantons dis- 
orderedjas it is vital for the considered above construction 
of Ref. §. 

The models considered in this paper include the scalars 
and fermions in a special way resolving the above men- 
tioned problems. The simplest model is based upon the 
SU (2) gauge group. The models provide the desirable ef- 
fective interaction between instantons making their identi- 
cal orientation the most probable. This result gives a hope 
that the phase transition into the state with polarized in- 



stantons is possible. Moreover, the mean field approxi- 
mation confirms the existence of the phase with polarized 
instantons. It is separated from the non-polarized phase 
by the first-order phase transition. 

It is important that the models considered provide the 
interaction of instantons belonging to the su{2) gauge al- 
gebra, but no interaction between the antiinstantons of 
this algebra. As a result one can expect no condensation of 
antiinstantons belonging to this su{2) gauge algebra. The 
models can describe the reversed situation: for a given 
SU{2) gauge group they can provide the antiinstantons 
with the interaction while instantons would possess no in- 
teraction. Then in this SU{2) gauge group one can expect 
the antiinstantons to become polarized, while instantons 
remain non-polarized. This is exactly what is necessary 
for the discussed above 50(4) gauge group construction, 
where instantons must be polarized for one su(2) subal- 
gebra and antiinstantons - for the other. 

Notice that for SO {4) gauge group the necessary vac- 
uum state of polarized instantons and antiinstantons does 
not violate the parity conservation law in spite of the fact 
that there is a clear distinction between instantons and an- 
tiinstanton. The polarized instantons belong to the su{2) 
subalgebra 1, the antiinstantons to the su{2) subalgebra 
2. The inversion results in the transformation P: instan- 
tons — > antiinstantons, and antiinstantons —>■ instantons. 
This, however, does not produce the new vacuum state. 
Really, there is a freedom to arbitrarily choose the two 
su{2) subalgebras. Simultaneously with inversion one can 
change the names of the subalgebras, P: subalgebra 1 
subalgebra 2, subalgebra 2 subalgebra 1. Under this 
transformation the polarized instantons remain in the sub- 
algebra 1, the antiinstantons in the subalgebra 2. In the 
new vacuum state the orientations of instantons and anti- 
instantons differ from orientations in the initial vacuum. 
In order to restore the initial orientation it is sufficient to 
fulfill the global rotation in the isotopic space, which is 
allowed due to gauge invariance of the theory. The vac- 
uum state remains invariant under thus understood inver- 
sion. If all the other properties of the theory are invari- 
ant under the transformation: subalgebra 1 — > subalgebra 
2, subalgebra 2 subalgebra 1, then the theory satis- 
fies the parity conservation law. In this paper we deal 
mainly with SU{2) gauge group. For this gauge group 
we are looking for the possibility to construct the vacuum 
with polarized instantons (or antiinstantons). The single 
SU(2) group with polarized instantons obviously violates 
the parity conservation law, because inversion results in 
the transformation of the instantons into antiinstantons. 
We will keep in mind that the parity would be restored, if 
one combines two SU{2) groups, having polarized instan- 
tons in one of them and polarized antiinstantons in the 
other. 

Consider the main idea. The usual way to govern the 
properties of the gauge theory is provided by the scalar 
condensate. It is clear that this condensate itself can- 
not resolve the puzzle of the different behavior of instan- 
tons and antiinstantons. In order to do this let us re- 
member that the instantons and antiinstantons interact 
differently with right-hand and left-hand fermions. This 
property manifests itself most strongly for fermion zero- 
modes: the instantons produce the right-hand fermionic 
zero-modes, while the antiinstantons give the left-hand 
zero-modes Realizing that we can develop the fol- 

lowing construction. Let there exist some condensates of 



the scalars. Let the scalar-fermion interaction has the 
scalar-pseudoscalar vertex proportional to 1 — 75. Then 
the right-hand fermions interact with the scalar conden- 
sate, while the left-hand fermions do not. As a result the 
right-hand zero modes of fermions provide a connection 
between the scalars and instantons. It permits the scalars 
to influence upon the instantons in a specific manner. We 
will see that this influence results in the desirable inter- 
action between the instantons. No similar connection ap- 
pears between scalars and the antiinstantons because the 
left-hand fermions do not interact with the scalar conden- 
sates. 

This construction may be described in the usual terms 
as a one-fermion-loop correction to the gauge field ac- 
tion. Calculating it we are to consider the interaction 
of the fermions with the scalar condensate as well as 
with the gauge field created by several instantons. In 
order to find the necessary effect one must choose prop- 
erly the scalar condensates. The condensates create the 
field V applied to the fermions. It is clear that this field 
must not commute with the gauge field created by instan- 
tons, [7V, V] 0, where V is the covariant derivative 
7V = 7,,(a^ - iA'^^{x)T''), T" = t'^/2, a = 1, 2, 3 are the 
generators of the SU{2) gauge group, and A°(a;) is the 
gauge field created by several instantons. Otherwise no 
connection between the scalar condensates and instantons 
would appear. If one wishes to consider homogeneous field 
V, as is usual, then the only way to satisfy this condition 
is to suppose that the field V depends on the generators 
of the gauge group: V ~ TU{1 — 75). As a result we are 
to introduce into the problem some additional vector U 
to multiply the vector of generators T. The constant vec- 
tor U not only looks ugly but makes no good, as can be 
verified. We are to find the vector whose averaged value 
is zero, but the averaged values of its powers could play a 
role: < U" >= 0, < WU'' >= (l/3)U'^Sab, •■■ ■ The way 
to do this can only be provided by an additional symme- 
try. We are to consider some "additional" SU{2) group 
and identify its generators with the necessary vector U. 
There are two ways to introduce this additional symme- 
try. It may be considered either as a global symmetry 
or as a local gauge symmetry. For the first case one can 
visualize group as "a fiavor group" or "a group of gener- 
ations". This possibihty was first discussed in Ref. ||l^ . 
For the second case the gauge group considered becomes 
as big as SU{2) x SU{2). Both possibilities and shown to 
provide the system of instantons with the desired proper- 
ties. There is meanwhile a substantial difference between 
these two realizations of the model. For the global addi- 
tional group the gauge field acquire a mass. In contrast, 
the local realization of the additional symmetry permits 
to keep the gauge fields for some SU{2) gauge subgroup 
massless. The latter property is desirable because there 
must remain the massless gauge field if we wish to con- 
struct afterwards the massless gravitons with the help of 
this field. 

Note that there is a clear and interesting analogy be- 
tween the discussed models and the phenomenon of fer- 
romagnetism. The problem of the phase transition into 
a state with polarized instantons resembles the transi- 
tion into a ferromagnetic phase. The instantons play a 
role very similar to the role of magnetic impurities in the 
ferromagnetic case. The fermionic zero modes resemble 
the atomic outer electrons. In the pure gauge theory the 



fermionic zero-modes for several instantons remain degen- 
erate in accordance with the index theorem , . This 
resembles the situation when atomic outer electrons are 
well-localized on the atoms and their energy levels are de- 
generate. The interaction of right-hand zero-modes with 
the scalar condensates results in the splitting of the zero 
modes. Therefore these scalar condensates play the same 
role as the crystal field which results in the splitting of the 
atomic energy levels and creation of a conducting band. 
The splitting makes zero-modes look similar to the elec- 
trons in the conducting band. The problem of instanton 
interaction looks very close to the origin of the exchange 
integral describing interaction of the magnetic impurities 
in the ferromagnetic theory. The exchange integral be- 
tween impurities appears due to their interaction with the 
electrons of the conducting zone. Similarly the interac- 
tion of instantons with the right-hand zero-modes results 
in the effective interaction between instantons. If some 
other impurities do not interact with the electrons of the 
conducting zone then they do not play a role in the fer- 
romagnetic state. Similarly, the antiinstantons which do 
not interact with the right-hand fermions do not play a 
role in the considered problem. 

The precise results of this paper are the following. The 
fermionic determinant det (F) is calculated in the one-loop 
approximation when fermions interact with instantons as 
well as with the scalar condensate. If we define 

det (F) = exp(-S'F) , (1) 

then Sp may be considered as an effective action, i.e. a 
contribution of the fermions to the action for the gauge 
field. It is found that it describes the interaction of the 
instantons. For the simplest case of two well separated 
instantons this interaction is found to be 

^^ = ^^^-7, (2) 

where is the value of the scalar condensate, / is the 
coupling constant of the scalar-fermion interaction, m is 
the mass of fermions, pi, p2 are the radii of the instantons 
and r is their separation, 7 is the angle between the direc- 
tions of instanton orientation. It is assumed in (|^) that 
Pi:P2 ^ r. Formulas (|^), (Q) describe the interaction 
between the instantons. The action (j^), first reported in 
Ref. ijl^ , has a desired minimum at 7 = for the identical 
orientation of instantons. Note that usually the radiative 
corrections renormalize physical quantities. In the case 
considered they provide the new phenomenon: the inter- 
action between the instantons. 

The mean field approximation for the ensemble of in- 
stantons interacting via action Eq.(H) is considered. This 
approach confirms the existence of the nontrivial vacuum 
state with polarized instantons. It reveals also the first- 
order phase transition which separates the polarized phase 
from the non-polarized. 

In Sections II the model based on SU{2) gauge the- 
ory with additional global SU{2) symmetry is considered. 
The interaction of the instantons is considered in detail in 
Sections III-V. In Section VI the mean field approxima- 
tion is developed. In Section VII the SU (2) x SU{2) gauge 
theory model is discussed. Throughout the paper we use 
the set of definitions presented in Section IX. The vital 
properties of the fermionic zero-modes in the field of sev- 
eral instantons and their role in the instanton-instanton 
interaction is discussed Sections X,XI. 



II. THE SU{2) MODEL 

Consider the SU{2) gauge theory. Suppose that there 
are two generations of fermionic fields in the fundamen- 
tal representation of this gauge group. The masses of the 
fermions are supposed to be equal and we will treat them 
as a doublet in the space of generations. Suppose also that 
there are scalars in the vector representation of this gauge 
group. Consider three generations of scalars, e.g. a scalar 
triplet in the space of generations. Let us introduce the in- 
teraction between the scalars and the right-hand fermions 
described by the Lagrangian 

C,f{x) = f^|J+{x)Mx)U\s^^^PB{x) . (3) 

Here / is a dimensionless constant of scalar-fermion inter- 
action, ipAix) is the fermion doublet, indexes A, B = 1,2 
label the variables in the space of generations, ^i{x),i — 
1,2,3 is the triplet of scalar fields. There is a freedom 
of choosing the matrixes W = U\g, i = 1,2,3 describ- 
ing the coupling between different generations of fermions 
and scalars. We choose these matrixes to be the triplet of 
generators of rotations in the space of generations 

W = U\b = '^ab/2, 1,2,3 . (4) 

Note that Euclidean formulation is used, see for example 
Ref. 

Notice that the Lagrangian Eq.(H) obviously violates 
parity conservation law. It can be restored if, alongside 
with the considered SU{2) gauge group, we consider the 
other SU{2) gauge group in which the Lagrangian has the 
form similar to Eq.(^) but with the positive sign in front 
of 75. 

The scalar fields ^i{x) are in the vector representation 
of the gauge SU{2) group, therefore 

<^,{x) = $,^a(2:)T'^ , (5) 

where T° = t°/2, a = 1, 2, 3 are the generators of SU{2) 
gauge transformations. Thus we have nine scalar fields 
$i^a(a;). Suppose now that their nonlinear self- interaction 
results in developing of the scalar condensate which has 
the following form 

{'^tA^))cond = 4'Sia , (6) 

where is a constant. One can always fulfill this condition 
choosing appropriately the self-interaction of the scalars. 
Then it follows from (^),(||) that there appears the field 
V, 

V = fcj^ ifu) 1^ , (7) 

which influence upon the right-hand fermions. 

Our goal is to calculate the fermionic determinant 
det {—i'yW — im — iV). It depends on the gauge field 
which stands in the covariant derivative 

^-f^W^^^^{d^-iAl{x)T'^) . 

We wish to evaluate the determinant when the gauge field 
is created by several instantons. The gauge field and the 
field V created by the scalar condensate do not commute 
[7V, 0, as follows from Eq.^. This makes the de- 
terminant to depend non-trivially on the field V The 



determinant is calculated below in the one-loop approxi- 
mation. This means that the gauge field and the scalar 
condensate field are considered as the given external fields. 
Let us present the fermionic determinant as 



det (— — im — iV) = 

det (-i7V - im) det (1 + GV) , 



(8) 

(9) 



where G denotes the propagator describing behavior of 
the massive fermions in the gauge field 



G = (7V 



(10) 



The first factor det (— «7V — im) in (|^) describes the 
known fermion behavior in the pure gauge field. It is not 
relevant to the effect considered. The important informa- 
tion is contained in the second factor 



det (F) = det (1 + GV) = exp (-Sp) 



(11) 



which we will discuss in detail. The determinant det (F) 
may be considered as a contribution of the fermions to the 
action for the gauge field Sp- 

The instantons are known to create the zero-modes of 
the fermions. These zero-modes play a crucial role in the 
following calculations. Therefore it is useful to distinguish 
them in the fermion propagator. With this purpose let us 
introduce the projection operator P onto the states of zero 
modes. It satisfies the conditions 



= P , (7V)P 0, Sp(F) = 2k 



(12) 



Here k is the number of instantons. The number of zero- 
modes in the considered case of two fermionic generations 
in the fundamental representation of the SU (2) group is 
2k. The propagator (nO) may be presented as 



G — Go + Gi , 
Go = P/m, 

Gi = (l-P)(7V + m)-i(l-P). 



(13) 
(14) 
(15) 



To simplify calculations let us consider the instantons so 
small that the condition 



mp <C 1 , 



(16) 



where p is the instanton radius is fulfilled. Then the 
fermionic mass m may be considered as a small param- 
eter and we will put m = wherever it is possible. Notice 
that the instanton separation is not restricted in this con- 
sideration. We will return to this point later, see Section 
V.A. 

The propagator of the nonzero-modes for massless case 
is simplified to be 



g(m-O) =(l_p)(^V)-l(l-P). 



(17) 



This propagator may be presented in the convenient form 
found in Ref. @] 



Gr°'-(7V)^^ + ^(7V) 
It is clear from (@),(|l8l) that 

yGi'"=°V = . 



1 1 + 75 , 1 



1 - 75 



(18) 



(19) 



Using Eq.(|T|) we find that only the zero-modes ( p^ 
give a contribution to (^ij) , while nonzero- modes are elim- 
inated 



det (F) = det ( 1 + —PV 
' m 



(20) 



This result greatly simplifies the problem because for the 
finite number of the instantons the operator PV is pre- 
sented by a finite matrix. 

It is necessary to keep in mind that the scalar conden- 
sate (^ results in the creation of the mass for the gauge 
field through the Higgs mechanism n& 



(21) 



where g is the gauge coupling constant. We will neglect 
this mass supposing that the instantons radii p are suffi- 
ciently small 



pMv < 1 



III. THE FERMIONIC DETERMINANT 



(22) 



In order to calculate the determinant ( |20D we are to 
consider the matrix V in the subspace of zero-modes 
where it has a form V = {t, B\V\s,A), here the indexes 
s,t = l,---,k numerate the zero-modes in one genera- 
tion of instantons, while the indexes A,B — 1,2 label the 
generations of zero-modes. According to (jsj) the wave- 
functions of the zero-modes have the form 

= *s.A(a;, a,n) = -[i+(a;)]„,„'e„',aWA • (23) 

Here the index a — 1, 2 is a spinor index of the right-hand 
spinors, the index n = 1,2 is the isospinor index, is 
the usual 2x2 antisymmetric tensor. Pt{x) is a quater- 
nion conjugate to Ls{x), defined in Eq. (|79|) . According to 
this definition [L+{x)]n,n' = Ls,p.ix){T~)nM' ■ The func- 
tion UJA in ( P^ ) describes the orientation of a zero-mode 
in the two-dimensional space of fermionic generations, it 
is normalized as (uj^Iuia) = Sab- 

Using the wave- functions (|2^) we find the matrix of the 
operator V (^ 

{t,B\V\s,A)^ J d^x^ls{x,a,n'){n',B\V\n,A) 

^'s,A(a;,a,n) , (24) 

where 



{n,B\V\n, A) 



(25) 



After sim ple algebraic transformations we evaluate from 
(|23|) , (p4[) , (|25|) the convenient form for the matrix 



where 



{t,B\V\s,A) ^-ti^Wt,sUBA , 



Wt,s = ^Re J Lt{x)iTLtix) 



(26) 



(27) 



s,t = 1, - ■ ■ ,k. According to (| 
symmetric 



the matrix Ws t is anti- 



Wt. 



-Wi 



t,s 



(28) 



Using (|26| ) we can find the exphcit expression for the de- 
terminant ( |20| ) for a small number k of instantons. For a 
trivial case of single instanton k = 1 we find det (F) — 1, 
due to Eq.(p8|). For two instantons, k — 2, we find from 
Eq.® 



PV = 



2 





«Tyi.2C/ 



6 



(29) 



Remember that U is the 2x2 matrix in the space of 
generations. In Eq.(p9|) is the zero 2x2 matrix in the 
same space. Diagonalize the matrix in Eq.(p9[) we find the 
determinant Eq.(po|) 



det (F) 



Parameter ^ is 



Zm 



(30) 



(31) 



It is instructive to look at Eq. (]30|) from the point of view 
of the eigenvalue problem for the operator 7V + V. In the 
absence of the perturbation there are four zero modes, 
which are the solutions of = with zero eigenvalue. 
The perturbation V results in the shift of the eigenvalues 
from the zero value 



(7V + V)iIj ^£1p . 



(32) 



We restrict our consideration to the subspace of the zero 
modes, where the operator 7V + is identical to V, 7V + 
V = PV, given in Eq.(p9|). The four- fold degenerate zero 
eigenvalue is split into two two-fold degenerate eigenvalues 

6+ and e_. 



£±=± — \Wi,2\ 



(33) 



Using them one finds the determinant det(F) (1 — 
£+/to)^(1— e_/m)^ which, of cause, is identical to Eq. (pO|) . 
This consideration clarifies the physical picture discussed 
in the introduction: the considered fermionic determinant 
appears due to splitting of the zero modes. 

An important case provides the perturbation theory 
limit C^W^s < 1- We find from (|ol),(|ll|) 



1,2 



(34) 



The presented results show that the fermionic deter- 
minant depends on the matrix elements Wt.s given in 
Eq. (p7[) . If all instantons have identical orientation qt = 
ptw, t = 1, • • • , fc, where pt is radius of i-th instanton, and 



IV. LARGE SEPARATION OF TWO INSTANTONS 

Consider the situation when the radii of the two instan- 
tons are less then their separation, and therefore the dilute 
gas approximation is valid. Let us calculate explicitly the 
fermionic determinant in this case. According to Eq.(p6|) 
we are to calculate the matrix element Wi^2 defined in 
Eq.(p7|). The functions Ls{x) in the dilute gas approxi- 
mation are given in ( pO| ) . Substituting them in Eq. ( p7| ) we 
find 



W"i.2 = ^Re 



yi - 2/2 ^4 

— —qiiTq2-, —d X 



yi 



2/2 1 



(35) 



The main contribution to the integral in the above expres- 
sion comes from the region \x — yi \ ~ |a; — ?/2| ~ I2/2 — 2/i| 
which justifies inequality (89) used to evaluate Lsjx ) (pO| ). 
The simple form obtained for the matrix element (|35|) is 
due to the chosen singular gauge which makes the vector 
potential to be strongly localized on the instantons, see 
Eq.(^). For our purposes it is sufficient to take the wave 
functions approximated by (^) in spite of the fact that 
they satisfy the orthogonal condition ( ^6| ) with accuracy 
^ l/r^2j which is comparable with the right-hand side of 
(|35|). This disadvantage does not manifest itself because 
we calculate in ( p5| ) the matrix element of the "imaginary" 
operator ir while in the orthogonal condition (|8^ ) the ma- 
trix element from "real" unity operator is calculated, see 
numerical results presented below in Fig.l. 
Calculating the integral in (^) we find 



1 



M^i,2 = — {Sp.-^ - '2Vf,v^)Re{T+q{ iTq2T^ ) 



(36) 



' 12 



where ru = \y2 - yi\ and = (1/2 - 2/i)/l2/2 - yi\- Taking 
into account identities r+qr" — 4:Re{q) and Re{vqv'^) = 



Re(g) valid for arbitrary quaternions g, v we find 



Wi,2 = —Re{q^iTq2) , 



W^f.2 = — ki|'l'?2|'sin2 7i2 , 
^12 



(37) 
(38) 



where 712 is the angle beween the directions of the ori- 
entation of instantons 1 and 2 defined by the identity 
Ke{qfq2) — |(?i||(72| COS712. Remember that this angle is 
defined up to mod tt because the instantons described by 
quaternions q and --q give the same field, see Eq.(66). 

Substituting Eq.(|3^) into Eq.(|30|) we find the determi- 
nant for two well separated instantons 



det(F) 



M'\<12\\ .2., 

C 4 sm 712 



(39) 



We are mainly interested in the region of large separation 



w satisfying condition w"* 



1 describes the instantons between instantons. If it is so large that C W^t s ^ li then 



orientation, then the matrix elements vanish Wt^s — 0. 
This follows from representation ( [79| ) for the functions 
Ls{x) given in Section VIII, but can be easily understood 
without calculations. For identical orientation there is no 
"complex" parameter in the problem. This eliminates the 
real part of the matrix element from the "complex" opera- 
tor if in Eq. (|^ . As a result the determinant in this case 
is trivial det(i^) = 1. We will see below that in general the 
determinant is the function of the instanton orientations. 



we can consider the perturbation theory limit (|34|). Using 
(|38| ) we find the effective action describing the interaction 
of 2 instantons 



Q _ cO o^2 lgir|g2| 



■ sm 712 



(40) 



This formula reproduces Eq.(^) proclaimed in the in- 
troduction. The generalization of this result to the 



case of 



straightforward, Sp 



several instantons is 

Note the pecuharity of det(F) (p9|). It possesses a node 
for specific values of parameters C, gi, (72, ?'i2, 712- If C > 1 
then this node may be situated in the region of applicabil- 
ity of the formula. The node has a clear physical meaning. 
When it happens then there is a "compulsory" creation of 
the four fcrmion-antifermion pairs produced by the two 
considered instantons. To sec this consider the instanton 
contribution to the amplitude of the pair creation. It has a 
multiple node coming from the determinant and the poles 
produced by the Green functions describing the propaga- 
tion of the created fermion pairs in the field of instantons. 
The node and the poles compensate for each other result- 
ing in the finite amplitude for the physical process. This 
effect is similar to creation of a pair of massless quarks 
in the field of an instanton in QCD, see for example |p^ . 
The difference is that in the considered case the fermions 
possess mass m > 0. Note also that the fermionic deter- 
minant is finite in the limit to = 0. The usual node of the 
determinant det(— — im) in the pure gauge field at 
TO = 0, see E g . m } , is compensated for by the pole arising 
from C in Eq.(|39[). 

Let us summarize the conditions under which the ob- 
tained results are valid. The radii p and separations of 
instantons r must obviously satisfy condition p < r which 
justifies the dilute gas approximation. Eq.(|l6|), mp < 1, 
guarantees the validity of the approximation based on con- 
sideration of zero-modes only. Eq.(40) was evaluated from 
Eq.(^9[) in the perturbation theory limit Sp < 1 valid if 
the separation is large enough (C|sin7|)^/^p < r. The 
radii of instantons must be small enough to protect them 
from being cracked by the scalar condensate. This means 
according to Eq.(p2|) that p < \/{g(j)). 

According to Eq.(B3) there is the interaction between 
the instantons. Eq. (ffO|) shows that it makes their identical 
orientation 712 = be the most probable. We did not use 
perturbation theory over the parameter C, = f4>/{2m). In 
this sense the problem of interaction between instantons 
is solved exactly. 



V. THE SUPPRESSION OF INSTANTON 
INTERACTION IN THE ASYMPTOTIC REGIONS 

It follows from from Eq. (^0|) that the effective action for 
two instantons depends on their separation as Sp ~ l/r"*. 
It is interesting to consider the action averaged over the 
positions of instantons. This averaged action has the log- 
arithmic divergence (S'_F)av ~ / d'^r/r'^ for small separa- 
tion r = as well as for large separation r = 00. The 
purpose of this section is to study what is happening in 
the vicinity of r ~ and r ~ 00. It is shown that the 
interaction of instantons is suppressed in both these re- 
gions. The origin of the suppression is different. For small 
separation the two-instanton interaction itself is strongly 
suppressed when r < pi+ p2. In contrast, there is no sup- 
pression in the two-instanton problem for large separation 
where Eq.(^o|) gives a correct asymptotic result for the 
two-instanton problem. The existence of other instantons 
changes this result drastically: the two-instanton interac- 
tion becomes suppressed exponentially in the presence of 
the other instantons. 



A. Small separation of two instantons 

Consider first interaction of the two instantons in the 
region of small separation ri2 ^ \qi\,\q2\. Simple analyt- 
ical results may be obtained in this region if we assume 
that 



^?2 » 



qiq2 sm 712 1 



(41) 



This inequality covers the interesting for us now case of 
small separation rj2 < I9192I provided the orientation of 
the two instantons is close enough I712I ^ 1. Remem- 
ber that we are interesting in the "polarized" instantons, 
therefore their identical orientation is of particular inter- 
est. Inequality covers as well the large separation 
region, which was considered above. It is shown in Sec- 



tion XI, sec Eqs.(117), that the matrix element describing 
interaction of the two instantons satisfying Eq. is 



1^1,2 = 2F{m,K2 



Re{q^iTq2) 



(42) 



' 12 



Eq. ([42|) differs from Eq.(|37|) by the factor F{ki, K2), where 
Ki = 2l(jj|/r, i = 1,2, and F{ki, K2) is defined in Eq.( |ll8| ). 

It is convenient to present the action Eq.(^) in the 
following form 



r' 



2|« |2 

92 .2 
4 sm 712 

12 



(43) 



Here K is the factor which shows how strongly the action 
Sp deviates from the asymptotic expression Sp defined in 
Eq.(|^). It follows from Eqs.(||),(|4|) that for the case un- 
der consideration this factor is K — F^(ki, K2). Its behav- 
ior as a function of ri2/{\qi\ + \q2\) is illustrated in Fig.l. 
For large separation ru » + |(?2| the factor K is close 
to unity, as is should be in order to reproduce the large 
separation result Eq.(^0[). The interesting thing happens 
for small separation ri2 — s- 0. H ere K is decreasing very 
fast. It is easy to find from Eq.(llS) the asymptotic be- 



havior F(ki, K2) rj2 (In l/r'12-l- const) , ri2 — > 0. Fig.l 
shows that this drastic decrease happens starting from 
fi2 < I91I + 1 92 1- Thus the interaction of the two instan- 
tons is strongly suppressed for small separation. Accord- 
ing to Eq. ( p3| ) the averaged action in the small-separation 
region behaves as (S'F)av J Kdr/r. The later integral 
is convergent, the main contribution to it gives the region 
in which ri2 > l^il -l- |92|- The other interesting feature 
illustrated by Fig.l is the fact that F^(ki,K2) is almost 
independent on the ratio of the instanton radiuses 1(711/1921 
thus depending mainly on ri2/(|(?i| + \q2\)- 

Up to now we have discussed the region covered by 
Eq. ( 4l|) , which gives the strong restriction on 712 for small 
separation of instantons. This was done in order to de- 
velop analytical calculations as far as it is possible. The 
general case of small separation and arbitrary orienta- 
tion of two instantons can be handled by direct numer- 
ical calculations. The useful representation for the ma- 
trix element W'i_2 is given in Section XI in Eq.(llO). In 
general case the factor K defined by Eq.(^3|) is a func- 
tion of three variables K = -ftr(Ki, K2, 712)- The results 



of numerical calculations for K based on Eqs.(llC),([43|) 
are also presented in Fig.l. The suppression of configu- 
rations with the small separation ri2 < 1(711 + \q2\ is the 
general feature, it is valid for any orientations and radii 
of instantons. Notice the very smooth dependence of K 



on 7i2. It manifests itself mainly in the region of small 
separation, which is suppressed. Thus the considered fac- 
tor K with reasonably good accuracy may be considered 
as a function of only one variable ?'i2/(|9i| + |Q'2|), K Ki 
((|9i| + \q2\)/ri2, (|<7i| + M)/?'i2), where the function 
F {ki, K2) is defined in Eq.( |ll8| ). 

We conclude: the interaction of the instantons is sup- 
pressed by the factor K when their separation is small. 



B. Large separation of two instantons in the presence 
of other instantons 



Consider two instantons 1 and 2 separated by large dis- 
tance ri2 3> |(7i|,|(72|- Our goal is to examine what is 
happening with their interaction if there are other instan- 
tons in the vicinity of either instanton 1 or 2. We know 
that the most important role in the problem play the 
fermionic zero-modes. We assumed in the previous con- 
sideration that the scalar field is weak enough, it does not 
seriously change the behavior of the zero-modes. The ba- 
sis for this assumption provides the inequality fcf) <C 
When one considers very large separation of instantons, 
then the zero-mode is to be considered far outside of the 
instanton radius where even the weak scalar field can dras- 
tically change the zero-mode behavior. Note that in the 
latter consideration we continue to call the zero-modes 
those functions which originate from zero-modes of the 
pure gauge theory. The scalar condensate influence upon 
them, shifting their eigenvalues from the zero value but 
this shift is supposed be small. 

Calculating the fermionic determinant Eq.(||) one is to 
face an eigenvalue problem Eq.(|3^), (7 V + V) ^ = e^. It 
is convenient for our purposes to re-wright it as 



G{x - x', ehM^'M^V^^' , (44) 



where — — iy4^(x)r''/2, and G{x,e) is the Green 

fimction whose kernel G{e) is 



G{e) = 



1 



jd + V 



79-K)(l+75)/2 + £ 
92 + e(Vo - e) 



(45) 



where Vq = fcjyTU. Eq.(p4|) is convenient for the asymp- 
totic expansion because in the singular gauge the vector- 
potential A(x') is localized on the instanton. Consider 
the zero mode ipix) — ipsix) localized on the instanton 
s. Then the integration in Eq.(Q) is also localized in 
the vicinity of the instanton s. As a result we find from 
Eq.(E3) the asymptotic 



ipsix) ^ G{x - ys,e)Xs , \x - ys\ 



00 



(46) 



where = — / 7pA^(a;')'0s(x')(i''x'. It is easy to find the 
constant Xg. Really, we assume that the scalar field f(f> 
is so weak that the zero-mode is not infiuenced by the 
condensate in the instanton vicinity. Therefore may 
be found using the pure zero-modes, non-affected by the 
scalar condensate. The pure zero-mode (x) satisfies the 
asymptotic condition similar to Eq.(^) with the Green 
function being Go{x) = 'yd/d'^ = 'yx/{2TT'^x'^). Comparing 
this with Eqs.(p4[),(|90|) one finds that A^ is the left-hand 
Dirac spinor A^ = X^{a,n) = 2TTqs,nn'£n'ai^A ■ 



We see from Eq. (^5|) that the asymptotic behavior very 
strongly depends on the eigenvalue e. For the pure zero- 
mode e = there is only the possibility for the power- 
type decrease of the zero-mode with separation. Even the 
small shift of the eigenvalue from the zero value changes 
the situation drastically. If £ ^ 0, then there appears 
the constant s{Vo — e) in the denominator in the Green 
function Eq.([45|) resulting in the exponential decrease of 
the zero- mode Eq. (^) . 

This consideration justifies the results of Section IV for 
the two-instanton problem. Really, in this case there cer- 
tainly appear the nonzero eigenvalues, but their magni- 
tude decreases with separation as ~ l/r**. As a result 
one can neglect e « in the denominator of the Green 
function, thus resulting in the power type decrease of the 
zero-mode. 

The situation changes if there are other instantons. Let 
for example there is the instanton 1 in the vicinity of the 
instanton 1 and the instanton 2 in the vicinity of the in- 
stanton 2. Let their separations satisfy r^j^i,r22' ^ ''12- 
Then the eigenvalues of the zero-modes centered on the in- 
stantons 1, 1 are influenced mainly by their mutual inter- 
action, being independent on the large distance ri2. The 
same happens with eigenvalues of the zero- modes centered 
on the instantons 2, 2 . This results in the exponential de- 
crease of the zero-modes. The zero modes centered on 
instantons 1,1 or on 2, 2 decrease as 



ips ^ exp{~Ms\x - Xs\), s = l,2 



(47) 



Here = (e^ - EsVo)^/^. According to Eq.(|4|) it de- 
pends on the splitting Eg = Ss± which can be found from 
Eqs.(||),(|37|), e,± « ±/0|Q,||(j,,|sin7,,'/(2rf^')- de- 
pends as well on the eigenvalues fcfi/i, —3/0/4 of the op- 
erator Vq — f4>TU. The exponential decrease of the zero 
modes with separation results in the decrease of the ma- 
trix element W12 responsible for the interaction between 
instantons 1 and 2. This consideration shows that the in- 
teraction of the instantons 1, 2 can be described by the 
power law Eq. (Bol) only if 



ri2 < rn 



1 

70 



(48) 



We estimate here the splitting of the eigenvalues of instan- 
tons 1, 1 and 2,2 as \Ss± \ ~ f(j). 

We conclude, if ri2 > Tmax then the interaction between 
instantons decreases exponentially with separation ri2. 



VI. THE MEAN FIELD APPROXIMATION 

Let us examine the possibility of the phase transition 
into the state with polarized instantons in the mean field 
approximation. Consider the ensemble of instantons. Ac- 
cording to Eq.(^) their interaction is described by the 
action 



\qs\'\qt\' 



(1 - (?i.,^nt,^)^) 



(49) 



We took here into account that sin^ 7sf = 1 — {ns,fj,nt,fj,)^ , 
where n^^^ — qs,fj./\Qs\ is the unit vector characterizing 
the orientation of the instanton s. It is shown in Section 
V.A that the interaction of instantons is suppressed for 



small separation r^t < \Qs\ + \<lt\- In Section V.B it is 
shown that the interaction is suppressed for large separa- 
tion > ^ijiax as well, here r^nax is estimated in Eq.(p8[) 
''max ~ ^/{f'P)- The suppression of the interaction in the 
asymptotic regions is taken into account in Eq. ( ^ ) by ne- 
glecting the contribution of the suppressed regions. The 
prime over summation in Eq.(^) reminds one of this fact. 

The corresponding statistical sum for the ensemble of 
instantons is 



Z = w cxp(— 5f ) , 



(50) 



where summation runs over the number of instantons, 
their positions, radiuses and orientations. The summa- 
tion includes also the fluctuations above the pure instan- 
ton picture. The function w in Eq.(|50|) is the probability 
to find the noninteracting ensemble of instantons in some 
particular state. 

In the mean field approximation one instanton is to be 
considered in the effective field created by the other in- 
stantons. The corresponding action is to be obtained av- 
eraging the action ( ^ ) over the positions, radiuses and 
orientations of other instantons. Moreover, we are mainly 
interested in one and only one degree of freedom for the 
considered instanton, its orientation. Therefore it is natu- 
ral to fulfill the averaging over the radius of the considered 
instanton as well. As a result one finds the effective ac- 
tion in the mean-field approximation as the function of 
the instanton orientation n 



> , where A ■ 



2C^|r . (51) 



Here Pq =< > is the averaged square of the instanton 
radius. Remember that there exists the scalar condensate 
which makes this radius to be finite, po ^ l/My ^ 1/(50), 
see Eq.(|2l]). The quantity < n^Tii, > describes the aver- 
aged orientation of the instantons in the vacuum. For 
non-polarized vacuum, when there is no preferred orien- 
tation, < n^rii, >— 5^v/^- The factor describes the 
separation of the instanton from other ones summed over 
their positions 



1 



= N I (fr\ = 27r2iVln- 



(52) 



where N is the instanton concentration. We can put 
''mill ~ 2po ~ 1/(5'/'); because for shorter separation in- 
according to Eq. ( p8| ) 

n 

) « ln(.g//). Using 
Eq.(fc2|) one can present the parameter A in the form 



stantons do not interact, while rmax 
is rmax ~ l/(/0)- Thus ln(r,nax/''m 



A = SC'^vqNL = 8B , where B = (^vqNL . (53) 

Here vq = tt'^Pq/2 is the volume occupied by the instanton 
with the radius po and L = ln(rinax/''min) ~ ln(g//). 

The probability to find the instanton with orientation n 
in the mean-field approximation is w^{{n) = exp S'mf ('t-). 
The self-consistency condition for the averaged orientation 
< n^n^ > is 

f nnUi, exp(A < ncrTir > n^nr)dn . ^. 
< n^riy >= ^ — r ^ 4 ^ „ ^ ^ ^ ■ (54) 



/exp(^ < naUr > n„nr )dn 



The integration here covers the sphere n^n^ = 1. The fact 
of polarization of instantons means that there is at least 
one eigenvalue of the matrix < n^n^ > which exceeds 1/4. 



In order to find such an eigenvalue it is sufficient to use 
for this matrix the anzats 



2. , 



(55) 



in which it is supposed that the possible preferred orien 
tation takes place along the 4-axes. Substituting Eq.(p5| 
into Eq.(|5^) one finds the equation for 

z 



SB 



f{z), where z — 8Ba 



(56) 



and the function f{z) is 

4 J njexp{znl)dn 1 _ 1 Ii{z/2) - l2{z/2) 



3jexp{znl)dn 3 3 Io{z/2) - h {z/2) ' 

(57) 



Here Im{x), m = 0,1,2 are the modified Bessel func- 
tions. The function f{z) has the properties: f{z) 
1, z^oo, /'(O) = 1/12 « 0.083 < max(/(z)/z) « 0.110. 
The later inequality indicates that f"{z) changes the sigh. 
Fig. 2 shows f{z) as well as the left-hand side of Eq.(p6|) 
for different parameters B. The nontrivial solutions of 
Eq.(H) exist if 



B = C voNL >Bc^ 1.14 



(58) 



It is clear that vqN is to be small, vqN < 1, otherwise 
the instanton overlapping would diminish their interac- 
tion. It is natural to put aside the peculiar possibility that 
hi(g/ f) may be considered as a large parameter. Then the 
only way to satisfy inequality in Eq.(^8|) is to have C > 1- 
It is important that this region of C can be considered 
in the framework of the formalism presented above. Re- 
ally, the only place where the perturbation theory over 
C was used is the evaluation of the action Eq.(^0|). One 
should use here the perturbation theory over the param- 
eter C^p^pI sin^ 7i2/rf2 < 1, which can be small even for 
large C > 1- 

We see from Fig. 2 that if S > = 3/2 then the non- 
trivial solution is unique. In the region Be < B < B'^ there 
are two non-zero solutions. One of them is unstable in the 
following sense. The fluctuation of the instanton polariza- 
tion in the partition function described by the right-hand 
side of Eq. (p^) results in the stronger fluctuation of the 
averaged polarization in the left-hand side. The existence 
of nontrivial solutions means that there appears the state 
with polarized instantons. The fact that there are two so- 
lutions in the vicinity of the critical value Be indicates that 
the nontrivial phase appears as a result of the first-order 
phase transition. 



Notice that the action Eq.(49) depends on the scalar 



products of the vectors Ug^i^. This feature makes it look 
similar to the classical Heisenberg model for ferromagnet- 
ics. The only distinction is the second power of the scalar 
product {us^p.nt.^)'^ in Eq.(^). There is a clear topologi- 
cal reason for this peculiarity. An instanton describes the 
gauge field which belongs to the adjoint representation 
of the gauge SU{2) group. The different SU{2) matrixes 
±n^T^ have the same representatives in the adjoint rep- 
resentation. In this sense = — n^. Therefore all the 
quantities describing the instanton can depend only on 
even powers of n^. In the mean field approximation this 
second power {ns.i_i'nt^^)'^ strongly manifests itself. It re- 
sults in the change of sign of f"{z) which gives two non- 
trivial solutions for Be < B < B'^, and consequently the 
first order phase transition. 



Let us discuss shortly the most important assumptions 
which are made in this section. 

1. We consider the contribution of the instantons to the 
statistical sum. The quantum fluctuations in the vacuum 
above the instantons are divided into two parts. One of 
them is the one- loop contribution of the fermions. These 
particular fluctuations are of grate importance because 
they provide the interaction between instanton making 
their identical orientation be more probable. They are 
taken into account in Eq. (^) explicitly through the action 
Sp. Notice that the role of these fluctuations is enhanced 
by the parameter C > 1- AH the rest fluctuations present 
themselves in Eq.(pffl) implicitly in the function w. In the 
mean field approximation all these fluctuations manifest 
themselves through the averaged values vo,N,L. 

2. We assumed that the picture based on the dilute gas 
approximation for the instantons is valid. There are two 
reasons in favor of this assumption. First, there is the 
scalar condensate in the model which suppresses the in- 
stantons of large radius. Second, if there is an overlapping 
between some instantons violating the picture of the dilute 
gas approximation, then still we know that the interaction 
of these overlapping instantons is suppressed. Therefore 
only the configurations with non-overlapping instantons 
should play an important role. 

3. The instantons are taken into account explicitly, the 
antiinstantons are neglected. The reason in favor of 
this approach gives a fact that in the considered model 
the instantons interact and the antiinstantons do not. 
There is however an interesting problem: there exists 
the instanton-antiinstanton interaction . It is necessary 
therefore to examine the role of this interaction in the con- 
sidered vacuum. This problem remains outside the scope 
of the present paper. 

4. The interaction of instantons is considered as a sum of 
two-bodies interactions. The virial corrections caused by 
simultaneous interaction of three or more instantons are 
not considered. This approximation is consistent with the 
dilute gas approximation. 

5. The mean field approximation for the ensemble of in- 
stantons is implemented in the most simple way. Only 
one degree of freedom - the polarization of instantons is 
taken into account in the effective action. All other de- 
grees of freedom are averaged out of the action. 

In conclusion, the mean field approximation confirms 
the existence of the vacuum with polarized instantons. 
The transition to the polarized state in this approximation 
is the first-order phase transition. 



VII. THE MODELS WITH 5(7(2) x SU{2) C G 
GAUGE SYMMETRY 

Consider the gauge group G which is big enough to 
include the product of two SU{2) groups: Gi x G2 C G, 
where Gi = G2 = SU{2). Then the idea of the model 
developed in Section II may be presented in the other 
useful way. It is sufficient to consider in this case one 
generation of fermions which belongs to the fundamental 
representation of both Gi and G2. Let ^{x) be a scalar 
field in the vector representations of Gi and G2 



Lagrangian 



(59) 



1-75, 



C,f(x)^f^f+{xMx)^^^ix) . 
Suppose that the field <I>(x) develops the condensate 



(60) 



mx))cond = 0T1T2 



(61) 



Then there appears the field V which influence upon the 
right-hand fermions 



V = f<t>TiT, 



1-75 



(62) 



Consider now several instantons whose gauge field is trans- 
formed under Gi. Then we can identify vector Ti in 
Eq.(ll) with f in Eq.(0) and vector f2 with U. The fact 
that T2 generates the gauge group does not play a role 
here because the instantons under consideration belong 
to Gi, they are not transformed under G2. Therefore the 
fermion determinant for the considered case is identical 
to the determinant discussed in the previous sections, and 
Eqs.(|30|), (^4|),(|39|) remain valid for the considered model 
as well. As we know there appears the interaction be- 
tween instantons with tendency to make their orientation 
be identical, see (^). The same is true for the instantons 
belonging to G2. 

Up to this point there is the close similarity beween 
the considered model and the one discussed previously in 
Section II. But there is an important distinction. The 



scalar condensate (61) is invariant under the SU{2) gauge 
group generated by T = 71+72. Therefore the three 
vector fields transformed under this group acquire no mass 
from the scalar condensate 



Mv = 



(63) 



here 7i,T2 are the generators of Gi and G2. Let us in- 
troduce the scalar-fermion interaction described by the 



This is in contrast to the model of Section II where all the 
gauge fields possess the mass (pl|). 



VIII. CONCLUSION 

It is shown that the considered models provide the 
strong interaction between instantons making their iden- 
tical orientation to be more preferable. At the same time 
there is no such interaction between the antiinstantons. 
One can easily reverse the situation changing the sign in 
front of 75 in Eq.(^ or in Eq.(60). Then antiinstantons 
interact and instantons do not. The interaction appears 
from one-fermion-loop correction to the gauge action. 

We considered two models. One based upon the SU (2) 
gauge group. The other model is based upon larger gauge 
group Gi X G2 C G, Gi = G2 = SU{2). These two mod- 
els agree in supplying the instantons with the interaction. 
The advantage of the model based upon the wider group is 
the fact that it preserves the massless gauge fields. These 
fields are transformed under the gauge group SU(2) over- 
lapping with the groups Gi , G2 in which the interaction of 
instantons takes place. If the condensate of polarized in- 
stantons is developed in Gi and (or) G2 then we will enjoy 
the possibility to consider the interaction of the massless 
gauge fields with that condensate. This is exactly what is 
necessary for the construction of Ref . ||l[ . 

It is shown that the mean field approximation confirms 
the existence of the phase with polarized instantons which 



is separated from the non-polarized phase by the first- 
order phase transition. 

The obtained results give a hope that there exists the 
state with polarized instantons in gauge theory. 



results in the general self-dual gauge field for k instantons 



iF^,{x)Ty2 = [V^,V,], where = 



A'^{x). For single instanton 
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IX. APPENDIX A 

In this paper the following definitions for isospinor and 
spinor matrixes are used 



r± = (±zf,l), a± = (±Z(7,l). 
The t'Hooft symbols are 

The Dirac matrixes 7^ in the Eucledean space are 7^ = 
(— 27*^, 7(J'^), where 7^^ are the usual Dirac matrixes in 
Minkowsky space. The Eucledean matrixes satisfy the 
condition {7^,71/} — 25 ^i,. In the spinor representation 
they have the form 



7m 



< 




75 



-1 
1 



Consider the /c-instanton general solution for the 
revew see Let us introduce a quaternion as q = 
where = l,---,4 is an arbitrary vector. In this 

notation = ■ Consider the (fc -f 1) x fc matrix 
Ms^t{x), s = l,---,/c-|-l, t=l,---,fc which has quater- 
nionic matrix elements. M(x) is supposed to be a linear 
function of the quaternion of coordinates x = x^t^ 

M{x) = B-Cx , 

where B and C are independent 1) x fc quaternionic 
matrixes. They must be chosen so that the condition 



M+{x)M{x) = R{x) 



(64) 



be fulfilled for any x. Here R{x) is a non-degenerate kx k 
real matrix. The matrix C can be chosen to be 



Sst, s, t = 1, • • • , fc 



(65) 



Having M{x) one can find the fc + 1 quaternionic vector 
N{x) which satisfies the equations 

M+{x)N{x) = 0, N+{x)N{x) = 1 . 

Then the vector-potential defined in the quaternion rep- 
resentation as 



M{x) = f „ % ) , N{x) = ( ^ + ) u 

where Pq{x) = 1 + p^/ja: — yp. The quaternion q describes 
the instanton orientation, which may be considered as a 
unit four-dimensional vector n = q/\q\, and the instanton 
radius p — \q\, the position of the instanton is given by 
the quaternion y. The quaternion u satisfying condition 
jwp = 1 describes the freedom of the gauge transforma- 
tions. The condition u — 1 corresponds to the singular 
gauge |0| in which the vector potential has the form 



i ^ qr^q+ix ~ y)i, 



2po{x) 



\x - 



(66) 



Similar simple formulas are valid for several instantons if 
they are well-separated p -^r, where p is the radius of an 
instanton and r is a separation of this instanton from the 
others. Then the dilute gas approximation considered in 
Ref. is valid. In this approximation the first row of the 
matrix M{x) is given by k constant quaternions 



Mi,t(x) t^l,---,k 



(67) 



while the remaining matrix elements may be approxi- 
mated as 



(68) 



Here qt describes the radius and orientation of the t-th 
instanton, yt is the position of this instanton. The com- 
ponents of the vector N{x) are 



Ni{x) 



Ni+t{x) 



1 



1 



^Po{x) \x-yt\ 



x~yt + , . 

ql u, t = I, ■ 



, k, 



where 



Po 



(^) = 1 + E 



f-[ \x-yt? 



(69) 
(70) 

(71) 



The vector potential for k well separated instantons in the 
singular gauge u = 1 has the form 



M^) = - 



2pq{x) 



\x-ytV 



There is a particular case [Q, jlS), ||l^ which is inter- 
esting for our purposes: the identical orientation of the 
instantons 



qt = Ptw 



(73) 



Here pt is the real radius of i-th instanton and the quater- 
nion w, vj'^w = 1, describes the common orientation of all 
the instantons. It is known that in this case Eqs.(|67|) 
- ^2) are valid for any separation between instantons. 



X. APPENDIX B. THE FERMIONIC 
ZERO-MODES 

The purpose of this section is to evaluate the conve- 
nient form for fermionic zero-modes in the field of several 
instantons. The projection operator P onto the space of 
fermionic zero-modes was obtained in |^] in the following 
form 



P = ( 1-(7V);^(7V) 



1-75 



(74) 



The following simple arguments can be used to verify 
it. It is easy to check out directly that Eqs.([l^) = 
P, (7V)P — are fulfilled and therefore P gives a projec- 
tion into the space of zero-modes. In order to prove that 
it gives the projection on the hole space of zero-modes it 
is sufficient to show that Sp{P) is equal to the number of 
zero-modes. It follows from (^4|) that 



Sp(P) = Sp 



1-75 



- (7V)^ 



1 1 + 75 

V2 2 



and therefore 



Sp(P) = Sp(-75) . 



(75) 



Every zero-mode gives a unity contribution to the right- 
hand side. The nonzero-modes do not contribute to it. 
Indeed, for any nonzero-mode ip\{x) satisfying the equa- 
tion 



i(7V)V'A(a:) = \tl^\{x), 



(76) 



with |A| > the function tp^\{x) — '^5ip\{x) satisfies 
—i{j\7)'ilj-\{x) = —X^p-\{x) and is, therefore, orthogo- 
nal to 'tp\{x). Thus the right-hand side of Eq.(|75|) is equal 
to the number of zero modes. This concludes the verifi- 
cation of the fact that P defined in ( [T^ is the projection 
operator on the zero-modes space. 

The Green function — describing propagation of 
scalars in the field of several instantons was found in Refs. 
0, For the considered case of isospin-1/2 it reads 



D{x,y)=-^ = 



1 N+{x)N{y) 



47r2 



(77) 



Substituting this expression into Eq.(^4|), after lengthy 
but conventional calculations presented below, see 
Eqs.(^-(102), it is possible to present the kernel P(x, y) 
of the projection operator in the following simple form 

P(x,y) = ^L+(a;)C+ci^(l-af)i(y) . (78) 



Here L{x) is defined as 

L{x) = R-\x)C+N{x) 



(79) 



i(a;)T+, where k 



The quantities M{x), N(x), B,C, R{x) constitute the 
ADHM construction Ref. ||lO| defined in Section VII. Note 
that L(x) is the vector with quaternionic components 
L{x) = {Ls{x), s = 1, • • • , /c), Ls{x) = is,p( 
is the number of instantons. 

The kernel of the projection operator satisfies the equa- 
tion (7V)P(a;, y) = for arbitrary y and the validity of 
this equation does not depend on specific properties of 
L{y) in the representation (uq). Therefore we can obtain 



the expression for a zero-mode wave function by replacing 
L(y) with arbitrary x-independent spinor-isospinor x- H 
we define 



V' = cruV(l-'75)x 



(80) 



then we find from Eq. (|7q) the wave functions of the 
fermionic zero modes 



It follows from Eq.(pO|) that x- independent spinor- 
isospinor ip satisfies the following conditions 



j5'ip = —"01 '^'''4' — — SV' 



(82) 



Let us define a — 1,2 the spinor index of the right-hand 
spinor and n — 1,2 the isospinor index. Then Eqs.(p^ 
result in 



1 

T2' 



(83) 



where is the right-hand spinor-isospinor defined by this 
equation and e„^Q, is the usual antisymmetric tensor with 
ei,2 — 1- Using Eq.(|8^) one can present the zero-modes 
given in Eq.(pl|) in the more detailed form 



1 



^s{x) = ^six, a,n) = - \L+{x)\^.e 



(84) 



The spinor tp is normalized in Eq.(p3|) as {ip^\ip) = 1. 
This results in the usual normalization conditions for the 
functions 5*5(3;) 



(85) 



The easiest way to check out these normalization condi- 
tions provides the equation P^ — P for the projection 
operator. For the kernel P{x, y) ([78|) it reads 



{2ir^)-^L+{x)C+Calr+ ■ 

d^zL{z)L+{z)<j+T+C+CL{ym - 75)72] 



{2^^)-^L+{x)a+r+C^CL{y)[{l 75)/2] 
To satisfy this condition L(x) must obey 

(27r2)-i(c7>+) {t+ I L{z)L+{z)d' 
(C+C)-V+r+ . 



Using identity (cr+o-+)(T+gT+) = 4cr+T+Re(g), which is 
valid for any quaternion q wc find 



Z T, = 



-^Re J L+{x)Ls{x) d^x = (C+C);/ = S^t 



(86) 



Here we choose the matrix C to satisfy identity {C'^C)^^ = 
Sst, which can be fulfilled for general fc-instanton solution, 
see (|65|). The normalization ( p5[ ) follows from (|8^). 

Formula ( pT] ) is valid for one generation of fermions in 
the fundamental representation of SU{2) gauge group. 
Generalization to the case of several generations of 
fermions is strait-forward. 



The formula becomes much simpler if the dilute gas 
approximation is vahd. It follows from Eqs. (|64|) , (|67|) , (pS]) 
that in this case the matrix R{x) simplifies to be 



Rst{x) ^ Sst {\x - yt\^ + \qt\^) , s,t 



1, 



(87) 



Using Eqs.(|87|),(|65|),(|69|),([70|) one finds the simple repre- 
sentation for the vector L{x) (|7^) in the dilute gas ap- 
proximation 



Ls{x) 



Eqs.(p8D,(|81|) give the fermionic zero modes in the dilute 
gas approximation. For single instanton they are reduced 
to the known expression 

Though Eq.(^ was evaluated in the dilute gas ap- 
proximation, it proves to be valid for the other important 
case when all the instantons have identical orientation and 
their separations are arbitrary. This follows from the com- 
ment given after Eq.(73). 

Eq. ( ^8| ) may be further simplified if we are interested in 
the region far outside of the cores of instantons 



where it reads 



\x - y.s\ 



qtu, s = 1, • • • ,A: , 



4 'is 



(89) 



(90) 



providing the very simple expression for the fermionic zero 
modes (pl|). 

Let us verify representation Eq. ( |7^ ) . With this purpose 
consider Eq.([77[) from which one finds 



:V^N+{x)N{y) 



(91) 



{x - y)' 



This gives 

(7V)i?(x,y)(7V) = (47r2)-i 
[-^^^^N+{x)N{y) 



1 



(x - yY 



[x - y)' 



:V,,N+{x)N{y)\V, 



(v) 

Here the last operator acts on the y-coordinate — Vl . 
To simplify the formulae it is usufuU to apply this operator 
to the left 



hV)Dix, y)(7V) = (4^2)-i7^.7. • (92) 
29(^) {{x ~ y)J{x - yf)N+{x)N{y) + 

2((x - y),/{x - yf)N+{x){W^y^N+{y))+ 
d^J'\l/{x-yf)V^N+{x)N{y)~ 

{\/{x-yf)V^^^N+{x){V\y^N+{y)y 



The first term in the square brackets above yields 

5(x — y) as it follows from the simple identity —At:'^5{x — 

y) = It^lvdi^^'ii.x - y)f,/ix - yY] = jf,-fudf,d,,[l/ix ~ 
yY]- Then from ( p^ ) we find for the kernel P{x, y) of the 
operator P (|7|) 



P{x.,y) = [ATT\x~yY]-\{l-l5)l'Alt.l<^ 
(2/(x - yY){{x - y)^N+{x){V.N+{y)) + 
{x-yY{V^N+{x))N+{y)) - 
{V,N+{x)){V,N+{y)Y 



(93) 



Let us calculate the terms with derivatives in the right- 
hand side of (^). Note that for the considered self- 
dual field the potential may be presented as A^(a;) = 
N+{x)d^,N{x), see Appendix C. Therefore V^7V+(a;) = 
d^N+{x){l - N{x)N+{x)). Using the equality 1 - 
N{x)N+{x) = M{x)R-^{x)M+{x) one finds 



V^N+{x) = dt,N+{x)M{x)R'^{x)M+{x) = 
-N+{x)df,M{x)R-\x)M+{x) = 
N+{x)Ct+R-\x)M+{x) . 



(94) 



Using the identity M+{x) = M+{y) - (x - y)+C+ one 
finds from (Q 



V^N+{x)Niy) = ~N+ix)CR-\x)T+ix - y)+C+Niy) 

(95) 

because M^{y)N{y) ~ 0. From ( |95| ) one finds 

{x-y)^N+{x){V,N+[y))+ - (96) 
{x-y),{^,,N+{x))N+{y)^ 



N+{x)C 



{x - y)^{x - y), {R-\x) + R-\y)) 



iT^ [{x - y)^.vlxR~\v) ^{x- yYvlxR-\x)\ 



{x - y)> 



C+N{y) 



Here the identity r+r^ = + iv'^v'^'^ was used. Now 
we use Eq. (|96|) to simplify the first two terms in square 
brackets in ( p3|) and Eq.(Q) for the last term in (|9^). As 
a result we find 



P{x, y) = [47r2(x - yY]-\{l ~ l^)l^]l^.lu 



(97) 



L+{x) 



{l/2){R{x) + R{y))5^, 



2iT''{n^,nxT]lxR{x) - n^nxr^^^^Riy)) 



T+M+{x)M{y)T- 



L{y) 



Here is a unit vector n^ — {x — y)^/ \/ {x — yY- 

Now let us use the following identity for the Dirac ma- 
trixes 



7m7.(1 - 75) = ((^M- + «<.f^")(l - 75) 



(98) 



Substituting this identity into (|9^) we get two terms, one 
comes from the 5^^^ and the other one - from iJ/J^yCr" in the 
right-hand side of (98) 



P{x, y) = P(s) {x, y) -f P(„,,) [x, y) . (99) 

For the first term we find substituting 7^71^ — s- 5f^u in ( p7| ) 

P(5)(a;,y) = [47r2(x-y)2]-i[(l-75)/2]- (100) 
L+{x) [2{R{x) + R{y)) - t+ M+ {x)M{y)T;] L{y) . 



With the help of the identity t^qt^ = 4Re(g) which is 
vahd for any quaternion q one finds that 

T+M+{x)M{y)T- = 2[R{x) + R{y) - {x - yfC+C] . 



Then it follows from (|lOq ) 

P(S){x,y) = ^^—^L+{x)C+CL{y) . (101) 



The second term _in_ ( p9[ ) is to be found substituting 



li^lu iJ/Lcr'^ in 



P(u,a){x.y) = [4^2(x- y)2]-i[(l -75)/2] • 
L+{x)inl,<7''{~2iT\n^n,^rt^R{x) - n.nxrilM] " 
T+M+{x)M{y)T-}L{y) . 

The first term in curly brackets is evaluated with the help 
of the identity r]1^u'r/l\'n^n\ = ^5ab- The second term is 
simplified with the help of ij'j^iyT^qT^ — 4iT'^Re(g) where 
q is an arbitrary quaternion. As a result one finds 

Pi^r,a)ix,y) - - ^ V"L(y) . (102) 



Combining (^ , (|lOl|) , (p^ we justify Eq.(|7|) due to 
identity crT^rJ = 1 — af. 



the 



XI. APPENDIX C. THE MATRIX ELEMENT Wia 

In this section the explicit convenient expression for the 
matrix element Wi^2 Eq.(p7|) is presented. Introducing 
convenient notation W = if- Wi^2 one can present Eq.(^ 



W = -^Im j L+{x)L2{x)d^x = 



(103) 



^Im / {L+{x)L2{x) - L+{x)Li[x)) 



Here the imaginary part of a quaternion q is defined as 
Im q = {q — q'^)/2. The second equality in Eq.(103) fol- 
lows from the antisymmetry condition ( |28| ) of the matrix 
element W12 = —W2i- Using Eq.([79|) for the functions 

Li{x), i — 1,2 one finds 

W = -^Im J N+{x)CR-\x)^R'\x)C+N{x)d'^x , 

(104) 

where the matrix E is defined as 



1 

-1 



(105) 



For two instantons the matrixes M{x),C in the ADHM 
solution have an explicit simple form 



M{x) 



91 92 \ 








yi-x b 


c- 




• 


b y2-x 









(106) 



where qi, i = 1,2 are the quaternions describing the radius 
and orientation of the two instantons, yi are the quater- 
nions describing the positions of the instantons and 



2/12 

2|yi2| 



(9^91 - 9^92) , yi2 ^yi-y2 



(107) 



To simplify Eq.(104) it is convenient to express the vector 
N{x) in terms of the matrix M{x). With this purpose let 
us use the identity 

N{x) = gnNq, Il{x) = 1 - M{x)R-\x)M+{x), (108) 

where Tl{x) is the projection operator on the vector N{x), 
No is any vector non-orthogonal to N{x): N'^{x)No ^ 0, 
and UN is the normalization coefficient which satisfies 
|ajv|2 = {N+{1- M{x)R-'^{x)M+{x))Noy\ It is con- 
venient to chose 



N, 



-(1,0,0) 



(109) 



Substituting Eqjl0|) , (|To|) in Eq.(|To|), using the explicit 
representation ( |106D for the matrix M{x) and the obvious 
identities C+Nq = and R-^{x)Y,R-^{x) = i;/deti?(x) 
one finds 

1 f QR-\x)m+{x)j:m{x)R-^{x)Q+ . 
W = —Trim / ^tt; — , , , „ , , a X 



detRix) (1 - QR-^{x)Q+) 



(110) 



Here m{x) is "the square part" of M{x): 

and Q, are the two-dimensional quaternions 

Q = (91,92) , 



9j 



(112) 



Representation (IIC) is convenient for direct numerical 
calculations as well as for asymptotic expansion in differ- 
ent regions. The most simple and interesting is the case 
of large separation, — \yi2\ ^ |9i|, |92|- For this case 
the main contribution to the integral in Eq.(llO) gives the 
region \yi — x\ \y2 — x\ ^ ri2 in which 



m{x) « mo (a;) 



yi- X 
y2~ X 



and 



Rix) 



\yi - x\ 





12/2 



(113) 



(114) 



Substituting Eqs.(|ll3|),(pl4|) in Eq.(pIC|) and taking into 
account that \QR^^{x)Q^\ <C 1 one finds 



W = 2 



Im(gig^) 



(115) 



' 12 



which agrees with Eq.(p7|) evaluated in Section IV using 
the different technic. 

The other important case provides the region ^ 
1(7152 sin 71 2 1 in which the instanton separation might be 
not only large but small as well, less then the instantons 
radiuses, provided th eir o rientations are close I712I <Si 1. 
In this case b in Eq.(lll) is smaller then the separation 
|6| < C r i2. Therefore one can again use the asymptotic 
Eq.( |ll3| ) for ■m{x). In contrast the matrix R{x) is to be 
considered more accurately. 



R{x)tj w q+Qj + {m+ {x)moix)) ^ 



i,j = l,2, (116) 



because the separ a tion might be l ess then the radiuses. 
Substituting (|ll3| ),(116) in Eq.(llO) one finds after simple 
algebraic calculations 



(117) 



' 12 



where Ki = 2\qi\/r, i = 1,2 and the function F(ki,K2) is 

F(«;i,«2) = ^ J {\x\^-l) X (118) 

d^x 



IP x+ 1 ' 



+ kI) i\x + + 4) - kIk. 



2^,21 

2\ 



The explicit form of this integral representation for 
F{ki, K2) is evaluated in the convenient coordinate frame 
yi = —y2 = 1. The integral in Eq.(llS) is easily reduced 



to the one-dimensional one. Its properties are discussed 
in detail is Section V. 
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Fig.l 




FIG. 1. The interaction of instantons for small separation. The coefficient K{X) describes how strongly the interaction of in- 
stantons deviates from the asymptotic value, seeEq.(43), is shown vs the relative separation between instantons X — ri2/(pi+p2). 
The different curves present the dependence of K{X) on the ratio of instanton radiuses pi/p2 and on the angle 712 between 
the directions of the instanton orientations, cos 712 = n-i.^n2,^i. The full curve: pi/p'z = 1, 712 = 0, the dashed curve: 
pi/p2 = 5, 712 = 0, the dashed-dotted curve: pi/p2 = 10, 712 = 0, the dotted curve: pi/p2 = 1, 712 = 7r/2. The coefH- 
cient K{X) is found as K{X) = F'^{ki, K2), see Eq.(118), for 712 = 0. For calculations with 712 = 7r/2 the more general Eq.(llO) 
is used. The drop in K(X) for small separation indicates that in this region the instanton interaction is suppressed. 



Fig.2 




FIG. 2. The mean-field approximation. The full curve: the right-hand side of Eq.(56) defined in Eq.(57). The straight lines 
give the left-hand side of Eq.(56) for different parameters B. The dashed line: B < Be, the dotted line: Be < B < B'^, the 
dashed-dotted curve: B'^ < B, where Be = 1.14, B'^ = 3/2 and the three considered values of B are B = 0.75, 1.19, 1.75. The 
points "a,b,c" show the noiitrivial solutions in which the instantons arc polarized. The point "a" indicates the unique nontrivial 
solution exiting for large B. The points "b,c" are the two solutions existing for intermediate B. One of the solutions, the point 
"b" , is stable, the other solution, the point "c" , is non-stable to the fluctuations of the instanton polarization. The mean field 
approximation confirms the possibility of the state with polarized instantons. 



